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Abstract 

The fuzzy disc is a discretization of the algebra of functions on the two dimen- 
sional disc using finite matrices which preserves the action of the rotation group. 
We define a ip'^ scalar field theory on it and analyze numerically for three different 
limits for the rank of the matrix going to infinity. The numerical simulations reveal 
three different phases: uniform and disordered phases already the present in the 
commutative scalar field theory and a nonuniform ordered phase as a noncommut- 
ative effects. We have computed the transition curves between phases and their 
scaling. This is in agreement with studies on the fuzzy sphere, although the speed 
of convergence for the disc seems to be better. We have performed also three the 
limits for the theory in the cases of the theory going to the commutative plane or 
commutative disc. In this case the theory behaves differently, showing the intimate 
relationship between the nonuniform phase and noncommutative geometry. 



1 Introduction 



Field theory on noncommutative spaces (for a review see [Ij) are a interesting arena to 
study the properties of noncommutative geometry [2], [3l H]. We will not review here the 
motivations for a noncommutative geometry. Among the applications of noncommutative 
geometry are fuzzy spaces. The original fuzzy space is the fuzzy sphere but other 
spaces have been built (for a review see p]). Field theory on fuzzy spaces are a way 
to study field theory on noncommutative spaces on a finite setting. The fact that the 
approximation is based on matrices makes them ideal for numerical studies. The most 
studied case is the fuzzy sphere, for a review see [3, [8]. Most of the interest for these 
investigation is the presence of the limit for which the rank of the matrices grow, and at the 
same time the radius of the sphere increases, thus recovering in the limit a noncommutative 
plane. Here and in the following by noncommutative plane we mean the algebra generated 
by the noncommuting variables x and y with commutator 

[x,y] = ie (1.1) 

This sort of plane is usually described by a noncommutative * product among the fields, 
such as the Gronewold-Moyal product O [10] • Other product are possible, and a com- 
parison among them is given in [TT]. For our purposes it is however better to think of 
the "quantization" of the filed as a map from functions / to operators /. As we shall 
see the fuzzy disc approximation will render the matrices finite, and hence amenable to 
numerical simulations. 

We will study then a quantized ip^ scalar field theory approximating field with N x N 
matrices, and perform three different limits, giving rise to the noncommutative plane, 
the commutative plane and the commutative disc. We are interested in particular to the 
phase transitions fr2\ of the theory as we change the parameters of the action. The 
quantity of interest are susceptibility and specific heat, as well as other order parameters, 
which we will describe below. 

Is worthwhile compare the present method with the simulations of such theories on 
the lattice. In general the simulations of scalar theories on fuzzy spaces are slower than 
in their lattice counterparts since the fuzzy models due to the self-interaction term ip'^ are 
intrinsically non local and for higher power of the self-interacting the number of operations 
to calculate each Monte Carlo step SS grows even faster. However, we can expect some 
advantages of the fuzzy approach with the simulations of other field theories, in particular 
the fermionic theories or super symmetric models [H]. Beside another recent field of 
application the matrix models are the Yang-Mills actions [151 IIHl [U] which appear the 
concept of the emerging geometry. 

The paper is organized as follows. In the next section will be introduce the fuzzy 
disc and its properties. The real scalar field theory on the fuzzy disc will be presented in 
section [31 The simulation and the result are in Sect. [H where also the phase boundaries 
between the three phases present and their scaling properties are derived. Finally, the 
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results are summarized and discussed in the conclusions. 



2 The fuzzy disc 

The construction of the fuzzy disc [HI [191 ISHl 1211 122] is based on a quantization of the 
function on the plane, i.e. the association to a function of an operator on an infinite 
dimensional Hilbert space F. Consider the Hilbert space generated by the countable 
basis !'?/'„) with n positive or zero. Consider the operator a and its hermitian conjugate 
acting on the basis as 



alV'n) = y nO \'lljn-l) , 

a^|^„) = V(n + 1)^ |7/>„+i) , (2.1) 

With the commutation rule: 

[a,at]=^l. (2.2) 

One recognizes of course the usual creation and annihilation operator of single harmonic 
oscillator, with a slightly unusual normalization (in the following we will revert to the 
usual one). In terms of the variable of (11.11) . 

Consider a function on the plane in terms of the conjugate variables z, z with z = x+iy 
we define a quantization map which associates to a function an operator / according to 
the following rule: 

oo oo 

fiz,z) = J2 fl^'^""^^ ^ / = E /^"'^^"^a". (2.3) 

m,n=0 m,n=0 

Here and in the following we will not discuss the issues of convergence of the series, we 
assume the coefficients are of rapid decrease. A more rigorous construction, based on 
the coherent states of the Heisenberg group and on a Weyl-Wigner map can be found 
in Ell 



This quantization associates an operator to a function, the inverse map is 

/=(^l/>)- (2.4) 

with \z) the usual coherent states a|z) = z\z). The function associated to the operator is 
called the symbol of the operator. 

More generally we can consider operators written in a density matrix notation: 

oo 

/ = E /-n|V'm)(^„|. (2.5) 
m,n=0 

To this operator is associated the function: 

zz z^'z"^ 

f(z,z)=e~^ fnm , (2.6) 

m,n=0 * 
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Where the relation between the Taylor coefficients f[f^ and the fmn is 

, ,,,, 
Ilk- 2^ Jl-qk-q <Qq ' l^-'J 

q=0 ^' 

After quantization if the plane using the map fl2.3p and its inverse we define the 
subalgebras of finite (A^ + 1) x (A^ + 1) matrices considering the functions with a truncated 
expansion, more precisely taking in account only the expansion terms if both n and m are 
smaller than a given integer A^. Keeping 6 fixed the limit N oo the truncated algebra 
tends to the noncommutative plane. This is obtained with the help of the projection 
operator 

AT 

^r=Ei^")(^"i- (2.8) 

n=0 

To this operator corresponds the function 

n=0 n=0 n=0 

This sum can be expressed in terms of incomplete gamma function and gamma function 
obtaining a radial function: 

In the limit — )■ oo if 6* is fixed and nonzero the symbol f l2.10p converges pointwisely to the 
constant function Pg^\r,(p) = 1 and recovering the noncommutative plane. Otherwise if 
the product A^^ is fixed equals to a real constant the limit for A^ — )■ oo of the (12.101) 
become: 

1 r <R 
P^''\r,ip)^ 1/2 r = R 
r>R 



Mr)- (2.11) 



Where the convergence still pointwise. This limit converges to a step function in the radial 
coordinate r in other words the characteristic function of the disc on the plane. Given 
the algebra A generated by a and a^, a sequence of sub algebras Ag^^ can be defined by: 

= PnAPn. (2.12) 

The effect of this projection on a generic function is: 

N 

TL — Tfl 

m,n=0 ^ ' ' 

namely a truncation of the series expansion fl2.6p . The full algebra A is isomorphic to 
an algebra of operators however the previous relations shows that ^4^^'' is isomorphic to 
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Mjv+i the algebra of (A^+ 1) rank matrices, on each sub algebra A^f^^ the symbol f l2.10p is 
then the identity matrix in every Mjv+i- The name fuzzy disc derives from the cutoff; for a 
fixed N every function with the same first (A^ + 1)^ terms are mapped in the same matrix 
losing all the informations of the higher orders. In addiction the symbols of the fuzzyfied 
function are still defined outside the disc of radius = N9, but they are exponentially 
damped outside the disc. 

It s convenient to use a related way is to define the fuzzy disc using an association 
that maps directly a function defined on the disc of radius R to a. {N + 1)^ matrix. Is 
well know that in polar coordinates z = pe*'^ a base for /(p, -i?) are the Bessel functions 
of integer order which are the eigenf unctions of the Laplacian on the disc with Dirichlet 
boundary conditions. If / is square integrable with respect to the standard measure on 
the disc dx = pdpd^ it can be expanded in terms of Bessel functions: 



hOO OO 



/(P' ^) = E E fnr^''"'J\n\ (Vw) ■ (2-14) 

n=— OO k=l 

There is also a fuzzy version of the angle variable in terms of the phase operator [23]. 

It is possible to define a fuzzy Laplacian the eigenstates of this fuzzy Laplacian will 
be the analogous of the standard Bessel in the sense that the fuzzy Bessel are a base for 
Mjv+i and their symbols tends to the standard Bessel for — > oo. The first step is to 
define the derivations in using the projectors. In A we have: 



Therefore we define 



{z\ f,a^ \z) 



e 



1 



(TV) 



(TV) 



{z\ f,a \z) . 

prfpr,' 



pr k) • 



(2.15) 
(2.16) 

(2.17) 
(2.18) 



These operators satisfy the Leibniz rule and are a derivations on each A^^K In the same 
way for the definition of the Laplacian operator we start from the exact expression on A: 



V^f{z,z)=4dAf{z,z) 



'\A [a, [f,a)]] \z) . 



(2.19) 



We define in each A 



and then 



vvr(.,.)=^(.ipr 



p^fprra^ 

prfpr.a^ 



pr\z), 



p. 



(TV) 



(2.20) 
(2.21) 
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Can be proved that the eigenvalues of this Laplacian are the solution of this equation: 



N+l-n 

E( 

k=l 



\N+l-n+k 



(A^+ 1)!(A^+ 1 -n)! (X 



{N)' 
n,k 



k\{n + k)\{N + l-n + k) \ AN 



0. 



(2.22) 



There is a sequence of eigenvalues labeled by n, k where n is the order of the correspond- 
ent Bessel function and k indicates the k-th zero. This eigenvalues can be calculated 
numerically from fl2.22p and they appears to converge for — > oo to the spectrum of the 
standard Laplacian defined on a disc with Dirichlet homogeneous boundary conditions on 
the edge. The eigenstates of the standard Laplacian are: 



{N) 
n,k 



A 



(TV)- 
n,k 



N-n 

^/a\{a + n)\ 



a=0 



E 

s=0 



A 



(iV)- 
n,k 



slis + n)!(a — s) 



a+n I 



(2.23) 

Where the eigenstates with the negative n are hermitian conjug ate of In the 

fuzzy approximation the continuum eigenf unctions e*"'^J|„| {^yX\n\,k^") are represented by 
matrices M^v+i but the max possible n are fixed by the dimension of the fuzzyfication 
n< N. 



3 Real scalar field theory on the fuzzy disc 

The previous construction allows us to define a fuzzyfied version of a field theory on a 
disc. For the real scalar case with Euclidean signature, with the (p'^ potential, the action 
is given by: 

S{^) = I (^VV + /i^' + ^V')d'z, (3.1) 

with V the Laplacian on the disc, n and A are respectively a mass and an interaction 
parameters. We choose this particular model because for the plane case it is well studied. 
It is known [25] 126] that the diagrammatic expansion of this theory has only one divergent 
diagram, the so called tadpole diagram and it is renormalizable. Using the fuzzyfication 
procedure the action (13.11) can be approximated by: 

Sn^ = 710 Ti + + ^(^^^ = TT^^Tr (^^^ [a, a^]\ + ^i^^ + ^^^^ . (3.2) 

In this action (p are finite dimensional hermitian matrices and the product between the 
field became the standard matrix multiplication, being a finite matrix model can be 
approached numerically using Monte Carlo techniques a method currently in use for other 
fuzzy spaces like the sphere. Note that the fuzzy approximation provides both an infrared 
and an ultraviolet regularization of the theory. 

We can rescale the operators a,a^ (12.11) as a = \/9di,a^ = -\/^a^ in order the extract 
the 9 dependence and recast the action making 9 manifest: 

Sm{v) = vrTr (4^ [a, [<^, a)]] + 9 ( + ^0'] ] . (3.3) 
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It will be useful later to define the potential part of the action given by (from now to 
simplify notations we will omit the hats and the subscript A^, as all fields are matrices): 

Sv{v)=7^0Tr(^f^^' + ^v'y (3.4) 

In the simulation we will use another equivalent choice of the parameters often used for 
the potential part namely: 

5U^) = 7rerTr(^/iV + ^^. (3.5) 
The two choices are related by 

r = A, 

r/i' = /i. (3.6) 

Since the scalar field on either commutative or fuzzy disc is defined on a finite volume it 
is impossible to have a phase transition, however phase transitions may be found when the 
matrix dimension or the radius of the sphere become infinite. As already states the fuzzy 
disc arise from the relation = ON where 9 is the parameter of non-commutativity and 
R is the radius of the disc. Introducing in the non-commutativity parameter a dependence 
on the matrix size and performing different limits we have three different cases: 

CD — 7- oo with R^ = N9 fixed. In this case we recover the commutative disc of radius 
R. In the following we will take R = 1. 

NCP — 7- oo with i? — 7- oo and 9 fixed. In this case we recover the noncommutative 
plane. We will take 9 = 1 and therefore R^ = N. 

CP A^ — )■ oo, 6' — )■ and R — )■ oo. In this case we recover the commutative plane. In 
the following we will the 9 = l/y/N which implies R^ = Vn, in the limit N ^ oo 
we recover the commutative plan^. 

In the simulation we expect at least one critical line corresponding to the critical 
behavior of the continuous scalar field on the plane. Indeed one of the main aim of the 
simulation will be to find the phase diagrams for the fuzzy field theory. 



4 Simulations 

Using Monte Carlo method we will produce a sequence of configurations {4'j}j=i,2,--- ,Tmc 
and evaluate the average of the observables over the set of configurations. These sequences 

^It is possible recover the commutative plane taking 9 = TV", where a is a real number varying between 
— 1 < a < 0, our choice is just dictated by simplicity. 
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of configurations, called Monte Carlo chain, are representatives of the configuration space 
at given parameters. In this framework the expectation value is approximated as 



-1 Tmc 

{0)^-^Y.O^^ (4.1) 

where Oj is the value of the observable O evaluated in the j-sampled configuration, ifjj, 
Oj = 0[ipj]. The internal energy is defined as 

EiN,ix,X) = {S) , (4.2) 

It is also very useful to compute separately the average values of the two contributions: 

D{N,fi,X) = {TT{Sn)) , (4.3) 

ViN,fi,\) = {TiiSv)) . (4.4) 

Where So is defined as: 

SD{v) = 47rTT{ip[a,[ip,a^]]}. (4.5) 
The specific heat takes the form: 

C{N,fj.,\) = {S')-{S)K (4.6) 

These quantities correspond to the usual definitions for energy 

1 87 

W/^,A) = --— (4.7) 

and specific heat 

C(iV,/i,A) = — , (4.8) 

where Z is the partition function. Another important observable used to look for a phase 
transition is finite-volume susceptibility defined as: 

X=(Tr(y,)2)-(|Tr(^)|)^ (4.9) 

The phase transitions are located at the values the previous quantities show a peak when 
varying the parameters. Note that in comparing the 04.91) with the usual susceptibility 
the second term is contains (| Tr((^)|)^, otherwise due to the fact that symmetry causes 
Tr(y9) = there will be no peak even in presence of a phase transition. 

The simulations were conducted following the standard way using a Metropolis al- 
gorithm with the jackknife and binning methods [271 EH] in order to evaluate the errors 
on the expectation values (14.11) . The usual update algorithm, used to generate a new 
configuration, changes all the entries of the independent matrices: 

iPij — > ipij + ttij, with i,j = l,---N. (4.10) 
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where N is the matrix size and aij are random real numbers. The amphtude of variation 
of each matrix entry is tuned dynamically during the run time, requiring an acceptance 
rat^ between 10% and 30%. The new proposed configurations are judged by the met- 
ropolis algorithm, which requires to compute the difference between the action for the 
new proposed configuration and the action for the configuration at the current time step 
6Sij{x) = S{ipij) — )■ (fij + x). The probability of accepting the new configuration ipij + x 
is calculated as min(e~'''^, 1). In order to reduce the computation time, instead a global 
matrix change, we have changed just one matrix entry for each Monte Carlo step. The 
coefficient to update is not chosen randomly but following an order used to vary all the 
coefficients after A^^ time steps. The interaction term proportional to (f"^ of the potential 
in the action (13. 3 p couples a matrix entry cpij to all the other elements on its line and 
column. This means that using the previous described update procedure, the computa- 
tion of the variation of the action SS (and consequently the Metropolis check) needs a 
number of operations which increases linearly with the matrix sizqj. The drawback of 
this approach are mainly two: the complication of the code and the dependence of the 
correlation time r in N. In fact two configurations will share at least one coefficient until 
A^^ Monte Carlo steps, thus this correlations introduce an increase of the correlation time 
increasing and presumably r is proportional to A^^. However, using this optimization 
and the implementation of parallel computing we were able to execute our simulation 
for A^ = 20 with a good precision and in reasonable time despite our limited computa- 
tion resource. The initial conditions of the the Markov chain are chosen randomly in 
the configuration space although we used two general types; hot initial conditions, which 
are configurations far from the minimum and cold start conditions, which correspond to 
configurations close to the minimum. As usual we have excluded all the time step before 
the thermalization process take place. 



4.1 Order parameters 

The previous introduced quantities are not sufficient if we want to measure the various 
contributions of different modes of the fields to the configuration (p. We need some control 
parameters, usually called order parameters. Recalling the fuzzyfication procedure we 
associate to a function on the disc a matrix given by coefficients fnm in the density 
matrix base expansion of (12. 5p : 

°° ^n^m 2 °° „n+m A'd(n—m) 

ip{z,Z)=e 9 y(pnm^==== = e y (pnm^====, (4.11) 
m,n=0 * m,n=0 * 



^The acceptance rate is defined as the number of new configuration accepted over the total number 
of configurations tried 

•^The computation of quadratic part of the action, which couples each element of the matrix with 
another one and the calculation of the observables using the formula (|4.1[) are always sub dominant in 
terms of the number of operations. 
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Due to the linearity of the expansion and since the action is even, the coefficients (fnm have 
the property to average to zero. As a first idea, we can think about a quantity related to 
the square norms of the fields, for example the sums X^nm Iv^nml^- This quantity is called 
the fuU-power-of-the-field; it can be computed as the trace of the square: 

= M^') (4.12) 

The quantity {ipa) alone is not a good order parameter because it does not distinguish 
contributions from the different modes, but we can use it as a reference to define the 
quantity: 

N 
n=0 

it is easy to see that the quantity (I4.13P is connected with the purely radial contribution, 
thus can be used to analyze the pure radial contribution to the fuU-power-of-the-field. We 
will see that (ipl) ~ defines a disordered phase while {(pi) > defines a kind of ordered 
regime. In order to describe the type of ordering it is necessary to study the contribution 
from higher modes to ip"^. 

We can generalize the previous quantity and define parameters (pi in such a way that 
they form a decomposition of ip^: 

l>0 

Following this prescription, the other quantities for / > can be defined as: 

N-l 
n=0 

If the contribution is dominated by the radial symmetric parameter we expect to have 
{'-Pa) ~ (v'o)- ■'■^ next simulations we will evaluate, apart from / = 0, the quantity with 
/ = 1 as representatives of those contributions where the rotational symmetry is broken. 
According to fl4.15p we have 

N-l 

"^1 = ^ {\^n,n+l? + \Vn+l,n?) ■ (4.16) 
n=0 

Using higher I in fl4.15p we could analyses the contributions of the remaining modes, but 
it turns out that the measurements of the first two modes are enough to characterize the 
behavior of the system. 

4.2 Phase transitions 

Before we start let us remind of some aspects of the 2-dimensional y?^ model on a fuzzy 
sphere. Many simulations have focused on scalar fields on fuzzy spheres and they lead 
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show the presence of new phases among the classical ones, characterized as nonuniform 
(or matrix or striped) ordered phase |71 |29l [301 EH [32|, 133] . 

Numerical studies show three different phases: a disordered phase and an ordered 
regime, which splits into phases of uniform and nonuniform order. In the nonuniform 
phases, which is not present in the commutative planar ip^ theory, the rotational symmetry 
is spontaneously broken. The three transition curves intersect at the so called triple point. 
Analyzing the limit N oo limit (which for the fuzzy sphere leads to the noncommutative 
plane) is been found that the three curves and the triple point collapse using the same 
scaling function of N. Therefore, the triple point, all three phases and in particular the 
new nonuniform phase, survive in the limit. This means that the naive scalar field action 
can not be used as an approximation of the scalar field theory on the commutative 2- 
sphere. Using a perturbative approach it is possible to prove that in the limit of 
the commutative sphere the two point function is non affected by the UV/IR mixing but 
again the fuzzyfication does not commute with the commutative limit. Taking a different 
limit to the noncommutative 2-plane the UV/IR mixing reappears. So we can see the 
appearance of the nonuniform phase as a finite version of the UV/IR mixing. 

These nonuniform phases are closely connected with spontaneous symmetry break- 
ing [36l [37], this appearance of the new phase seems in contradiction with the Coleman- 
Mermin- Wagner theorem which states that there can be no spontaneous symmetry break- 
ing of continuous symmetry on 2-dimensional commutative spaces. However the theorem 
is strongly based on the locality of interactions therefore is no direct generalization of the 
theorem for the noncommutative case, where the theory is nonlocal. In the work of Gub- 
ser and Sondhi [3H], where those new phases were originally conjectured, it is argued that 
this translation non-invariant phase is possible only in dimensions d > 3. Nevertheless, 
subsequent numerical simulation, in particular of Ambjorn and Catterall [39], show the 
existence of such a phase even in c? = 2. 

Now we are ready to discuss the result of the Monte Carlo simulations and the various 
phases which arise in our model. In the following, for the first phase transitions, we use 
parameters A or of (13.41) . Fig. [1] shows the susceptibility (14. 9 p plots for a fixed values of 
A = 1, varying the matrix size and for all the limits described. The location of a maximum 
in the susceptibility is used to identifies the critical value related to the phase transition 
for the continuous model. In Fig. [T] the maximum of the three susceptibility in the 
= 20 case can be easily extrapolated for /ic negative equals to —1.608 ± 0.006; —2.30 ± 
0.01; —1.56 ± 0.01 respectively for the noncommutative plane (NCP), commutative disc 
(CD), commutative plane (CP) limit. In general we find the critical parameter /i^ for 
negative values and for — /i ^ A. 

This critical behavior shows up in the specific heat density as well, in fact the specific 
heat can be used as an alternative in order to determine where the phase transition take 
place. Actually both quantities displays the same critical value but, comparing Fig. [Hand 
Fig. [21 we find a small difference in the critical value of fic, likely due to a finite volume 
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Figure 1: Plots of the susceptibility in the different hmits CD, NCP, and CP(e ^ I, N-'^, N-°'^ ), for 
A = 1 and iV = 5, 10, 15, 20. 

effect. In the following we adopt the the susceptibility criterion in order to find the 
transitions curves because it is easier to estimate the maximum of susceptibility respect 
to estimate the critical parameter through the specific heat step. Beside using the same 
statistic the errors are smaller for the susceptibility than in the case of the specific heat 
since in general more statistics is required for the specific heat to have the same precision. 
From the slopes of the specific heats and of the the susceptibility we can infer information 
about the order of the phase transition. In fact this phase transitions turn out to be of 
the second order for all the three cases. 

In order to gain some informations on the nature of the phases involved in the trans- 
itions we look at the order parameters defined in the previous section. In Fig. |3] are 
plotted the expectation values of the order parameters. Looking for values greater than 
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Figure 2: Plots of the specific lieat density in the different limits CD, NCP, and CP, for X = 1 and 

TV = 20 

the critical point ^ > ^c^^ find the (yjo) and (99^) much smaller than the full power of 
the field {^pD with a weak dependence on /x. 

The behavior suggest that all the modes contributes effectively to the full power of the 
field as random fiuctuations. This type of phase is usually named disordered phase due to 
the contribution of the complete composition of the field by random fiuctuations around 
the null constant average field. For smaller values of /i with respect to respect the critical 
values /i < yUc we find a phase characterized by (yj^) ^ 0, the pure radial contribution 
(V9q) is dominant while {ip\) ~ thus only one mode (the zero one) is mainly present in 
the field decompositions governing its behavior. In contrast with the previously described 
phase, where the full powers of the fields are basically a sums of random fiuctuations, 
this ones is named uniform phase. In this regime the our fields are essentially diagonal 
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Figure 3: Plot of the expectation values of order parameters in the different hmits for for A = 1 and 

N = 20 

matrices and the kinetic action part of the action by the regions around its minima at 
(p± = ±a/— 2yu/ A. These results are completely consistent with what one would expect 
for the usual plane so we can state that all the three limits possess the disordered-uniform 
phase transition observed for a (p'^ scalar field theory on the commutative plane. 

Repeating the previous analysis for different values of A and finding the corresponding 
critical values /ic we can plot the transitions curves disorder-uniform; the scaled plots are 
shown in Fig. |H All the graphs possess good scaling properties even for small matrix 
sizes without strong finite dimensional effects. Each plot has a different scaling coefficient 
for the noncommutative plane limit the lines collapse for a scaling like the square root 
of the matrix size, for the commutative disc limit the best collapse take place for N~^'^'^ 
and for the commutative plane for A^"°'^. Furthermore the phase transition line is well 
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Figure 4: Plot of scaled uniform-disordered transition curves in the different limits for N ~ 5, 10, 15, 20 

approximated by a linear function therefore performing a fit of the data we obtains the 
approximations : 

l^^^^iN, A) ^ -0.29AAr-o-5, n^^{N, A) ^ -0A9XN-''-''\ i^^^{N, A) ^ -OASXN"''-^ 

(4.17) 

The previous hnear regression has been done on the highest matrix size data set N = 20, 
which is our best approximation to the large N limits. We see that the two commutative 
case, CD and CP behave in a similar way among themselves, and are different from 
the noncommutative plane limit. Beside the uniform and disorder phases the simulations 
show the appearance of a new phase, in this case it is convenient to adopt the other choice 
of parameters (13.51) . We find this new phase transition for large r. Fig. [5] show an example 
of this phase transition for fixed r = 10^ and varying /i. Again the critical points can be 
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Figure 5: Plots of the susceptibility in the different limits 9 = l,N^^,N^^-^, for r = IQ^ and N = 
5,10,15,20. 

easy located using the maximum of the susceptibihty. As we can see the locations of the 
peaks for all the plots are pushed to the left increasing in particular if we compare the 
plot concerning the non-commutative limit and the other two, we observe a bigger shift of 
the locations of the critical points to the left of the two commutative limit plots respect 
the noncommutative one. We will see this difference clearly in the transition curves. 

Looking at Fig. Othe new phase appears characterized, for /i < /ic, by {ipD < (ipf) ^ 
{ipD and {ipD ^ 0, this is true for all the limits we have investigated, therefore we are 
dealing with the same phase. In particular we notice that in the regions fi < fic the full 
power of the fields becomes almost linear, after various fits on the data we obtain for all 
limits: 

^ ^ h'^rcf,) . (4.18) 
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Figure 6: Plot of the expectation values of order parameters in the different limits for for A = 1 and 

N = 20 

The factor is a normalization. We can generalize this result introducing another para- 
metrization analogous to ( 13 .Sp in witch we reintroduce A but just as fixed parameter: 

Sviv) = n9r Tr (^fx^^ + ^(^^^ . (4.19) 

In this case we have for the correspondent full power of the field (14.181) the fit: 

(4.20, 

We recast our model in this way in order to show the analogies with the fuzzy sphere 
|29] . for which equation (I4.20p is very similar. The behavior of the order parameters 
suggests that in the new phase the eigenvalues of the field are in a minimum of the 
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potential. Following [29] the f l4.20p hint that the neighborhood of the minima of the 
potential mainly contribute to the expectation values. The minima of the potential [29] 
are N disjoint orbits of the form 



On = \I^U^ (In © U\Ue ?7(iV) / [f/(n) X f/(iV - n)], (4.21) 

where n < y and 1„ are the unit matrices. In the case of a negligible kinetic term, 
the orbit with the largest phase space volume will dominated the field composition. In 
the previous case of the uniform phase, the minimum of the potential is clearly found 
minimizing the kinetic term as well, but in the new phase there is not a direct hint that 
the kinetic term might be negligible. To test this assertion we have studied the partial 
contributions of the kinetic part D and the potential part V to the total energy. As Fig. [7] 
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Figure 7: Plot of the partial contributions V and D of the internal energy E in the different limits for 
for r = 10'' and N = 21. 
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shows, our results confirm that for jj < fic the transition the average kinetic contribution 
is much smaller than the modulus of the potential. However, this behavior is more marked 
for more negative /x- values in the uniform order regime as well, a fact that does not allow 
to interpret the rising of the this new phase caused only from the potential dominance. 
This new phase was already found for the same scalar field theory on the fuzzy sphere 
and is usually referred as non-uniform phase or matrix phase [TJ |29l [30] . 
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Fi gure 8: Plot of scaled transition curves in the different limits for N = 5, 11, 15, 21. 

Fig. [H] shows the scaled non-uniform-disordered critical lines, each curves scale in the 
matrix size with a different scaling factor; for the noncommutative plane limit the curves 
collapse for a scaling like for A^"*^'^, for the commutative disc limit the best collapse take 
place for A^~^'^ scaling behavior and for the commutative plane for N~^'^. Actually the 
scalings of of the commutative limits are not true scalings, in the sense that the collapse 
takes place for r large enough, in fact looking closely at the data set Fig. Owe notice that 
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the point towards the origin does not collapse. We will discuss this issue later. 
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Figure 9: Plot of scaled transition curves close to origin for 9 = N^^,9 = N^^-^ and TV = 5, 11, 15, 21. 

The asymptotic behavior for large r can be estimated algebraically [29i| from the scaled 
susceptibility curves and the critical point for the pure potential model, in fact for = 31 
and r — )■ +00 

. , 17.5 

/i, r ^ . 4.22 

r 

Furthermore using the scaling factors the asymptotic phase transition lines we have the 
approximations: 

/x^^P(iV, r) ^ ;,CD(^^ ^) ^ ^,^^{N, r) ^ ^^N-'\ (4.23) 

This large r behavior has an asymptote in zero and fits the scaled critical lines well for 
N = 31. 

In Fig. [To] we compare the non scaled transition curves of the noncommutative limit 
case and the other commutative ones, in both cases the difference became bigger for large 
commutative transition curves, the difference is much more evident for the commutative 
disc. The reason of this discrepancy is the explicit 6 dependence of the potential term in 
the action (13. 3p . The noncommutative limit is characterized by ^ = 1 making the model 
very similar to the one obtained using the fuzzy sphere in which the potential part has 
no dependence on the noncommutative parameter. Switching to the other to cases we 
introduce in the pre factor a dependence on the matrices size, increasing for fixed r the 
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Figure 10: Plot of the comparisons of transition curves between the noncommutative plane hmit and 
commutative disc hmit (left) and between the noncommutative limit and commutative plane limit(right) 
for iV = 5, 11, 15, 21. 

effect is to lower the potential dominance and graphically to lower the transition curve. 
For each we will find the previous Vc shifted towards lower values and since we have for 
all the cases an finite asymptotic behavior, in the commutative limits we can lower any 
point of the transition curves just increasing N. 

5 Conclusions 

The scalar field action on the fuzzy disc, in analogy with equivalent simulations on the 
fuzzy sphere, shows the occurrence of a new nonuniform phase, peculiar of the noncom- 
mutative models. We found the existence of three behaviors in the model described by 
(13. 3p which arise for all the limit performed. We studied the scaling properties of the 
transition curves disordered-ordered obtaining a scaling coefficient in each case, the non- 
commutative case has the same coefficient as the simulations on the fuzzy sphere scaling 
like the square root of the matrix dimension, the other two cases show different scaling 
behavior. Nevertheless looking at the order parameters in all three cases we can recognize 
the two phases as the disordered an the ordered ones. The critical line for the non- 
uniform-disordered phase transition is also found and the noncommutative limit scales as 
N~^'^ and grows asymptotically for large yU parameters. In the other two commutative 
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limits it is possible to consider a good data collapse using our matrix size approximations 
only for fi large enough. The asymptotic behaviors of are well approximated for large r 
by a pure potential model. We have confirmed that even for this models that the non- 
uniform phase is connected to the dominance of the minimized potential over the kinetic 
term via the phase space volume of the largest orbit. The geometry for a scalar field 
theory, defined on a fuzzy space, is defined by the Laplacian operator which correspond 
to a particular kinetic term of the action, thus the new phase should be insensitive on 
the underlying geometry. In particular, this argument justifies the similar behavior of 
the noncommutative limit comparing to the results found for the fuzzy sphere studied in 
[23 ED El]- A substantial difference from the fuzzy sphere case is in the commutative 
limits, in fact the 6 dependence of the potential part in each two cases modify the data 
collapse of the transition curves disordered non-uniform making a good collapse possible 
only for large r respect N. In the fuzzy sphere models the usual scalar field action cannot 
be used as an approximation of the scalar field theory on the commutative sphere due 
the persistence of the non ordered phase in the continuous limit. The problem is patched 
with an additional term used to increase the weight of the kinetic part or equivalently to 
decrease the influence of the potential term in order to suppress the anomalous phase [35] . 
In our model, considering the commuting limits, the action of the 6 on the potential term 
act naturally reducing the ratio between the potential part and the kinetic one, so we can 
surmise that in the continuous limits the non-uniform phase does not exist. The ability to 
differentiate the commutative and noncommutative cases is one of the main advantages 
of the method described in this paper. 
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